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Abstract 

It is argued that cosmological models that feature a flow of energy from dark energy to dark 
matter may solve the coincidence problem of late acceleration (i.e., "why the energy densities of 
both components are of the same order precisely today?"). However, much refined and abundant 
observational data of the redshift evolution of the Hubble factor are needed to ascertain whether 
they can do the job. 
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I. INTRODUCTION 



As is well known, the popular vacuum-cold dark matter (ACDM) model fits rather 
well most observational data, however it suffers from two main -and lasting- shortcom- 
ings, namely: the low value of the vacuum energy (about 123 magnitude orders below the 
quantum field theory estimation) and the coincidence problem. The latter refers to the fact 
that despite the vacuum energy does not vary with time and dark matter energy depends 
on redshift as (1 + zY both energy densities are nowadays of the same order. This is why 



very often the vacuum component is replaced by an evolving dark energy field -see [l, 2] for 
recent reviews. Nevertheless, this does not satisfactorily solve the problem since the initial 
value of the dark energy ought to be fine tuned up to obtain the desired result; it only 



substitutes a fine tuning by another 



2|. 



An encouraging approach that is receiving increasing attention ^] rests on the assumption 
that both components do not evolve separately but interact with each other, i.e., 

p,n + ^Hpm = Q , and + 3H{p^ + P^) = -Q . (1) 

Here, Q stands for the interaction term, and the subscripts m and (p denote non-relativistic 
(pressureless) dark matter (DM) and dark energy (DE), respectively. Clearly, this represents 
a step forward because for Q > (i.e., for energy transfers from DE to DM) the ratio between 
both energies, r = pm/ P4>i evolves more slowly than in the ACDM model. Effectively, in the 
latter, r = —3Hr while in the interacting case its time evolution of governed by 



3Hr 



w + 



(2) 



where w = P,j,/ P4> < —1/3 denotes the equation of state parameter of DE. Equation ([2]) 
follows from ([1]) and Friedmann's equation for spatially flat universes 

= k\p^ + P^) (/t^^SvrG). (3) 

As we will see below, for reasonable and suitably chosen interaction terms, r either tends 
to a constant or varies more slowly than the scale factor, a(t), at late times. This certainly 
alleviates the coincidence problem. It is worthy of note that if Q were negative (signifying 
an energy transfer from DM to DE) the problem would only worsen, moreover at early 
times the energy density of DE would have been negative, and, finally, the second law 
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of thermodynamics could be violated 0, Except in chamaleon models 6|, baryons as 
well as radiation usually do not participate in the interaction. The former because of the 
tight constraints imposed by local gravity measurements [3, 8|, the latter because otherwise 
the photons would not longer follow geodesies which would affect precise measurements of 
deviations of radar signals grazing the sun. 

It is not the aim of this article to solve the coincidence problem -we believe it cannot be 
done at present- but to devise a simple strategy to do it, within general relativity, as soon 
as sufficiently accurate and abundant data on the evolution of the Hubble factor becomes 
available. This paper extends and generalizes a previous work of the authors {9]. 

Most studies of late acceleration implicitly assume that matter and dark energy interact 
gravitationally only. In the absence of an underlying symmetry that would suppress a 
matter-dark energy coupling (or interaction) there is no a priori reason to discard it. Further, 
the coupling is not only likely but inevitable 1^ and its introduction is less arbitrary than 



otherwise. In fact, the accretion of DE by black holes H 



12l | implies a non-vanishing and 



positive Q. While its small value modifies only very little the dependence of the energy 
densities on redshift (and cannot substantially alleviate the coincidence problem), it sets a 
minimum value for the interaction. 

The interaction may bring some further consequences. It can push the beginning of 
accelerated expansion era to higher redshifts. It can also erroneously suggest an equation 
of state for the dark energy of phantom type (i.e., w < —1), see ISj and references therein. 
Moreover, the interaction may give rise to fiuctuations in the count of galaxy clusters with 
redshift jl4|. Likewise, the coupling can alter the isothermal Maxwell-Boltzmann velocity 
distribution of weakly interacting massive particles in the galaxy halos whereby the 
average dark matter velocity is expected to increase significantly. Since the detection rates 
of the different experiments searching for dark matter strongly depend on this velocity it is 
crucial not to dismiss a priori the coupling [l^ . 

Ultimately, the existence or non-existence of the interaction is to be discerned observa- 
tionally. Whereas the current empirical data are insufficient to settle this issue, a couple of 
studies seem to favor the former possibility: {i) As it should be expected, the interaction 
alters the time required for a self-gravitating, collapsing, structure to reach equilibrium as 



well as the equilibrium configuration itself. Therefore the Layzer-Irvine equation 



17 



18| is 



to be generalized to take into account the interaction. In this regard, from the analysis of 
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the dynamics of 33 relaxed galaxy clusters -for which reliable x-ray, weak lensing and 



opti- 



cal data are available-, it has been reported a small but not vanishing interaction [19|, |20 |. 
(ii) Since the interaction modifies the evolution rate of the metric potentials the integrated 
Sachs- Wolfe component of the cosmic microwave background radiation (CMB) is enhanced. 
The cross-correlation of galaxy catalogs with CMB maps also hints at a small interaction 



2l|. 



This paper is organized as follows. Next section introduces three different expressions for 
Q and show that all three lead to a late, everlasting, era of either constant or quasi- const ant 
r regardless of whether w is lower or higher than —1. Section 3 shows the consistency of 
the three interacting models with the recent observational bounds of Daly et al. 2^ based 
on data from supernovae type la, taken from different sources, and 30 radio-galaxies. The 
only cosmological assumption behind these bounds is that the Universe can be considered 
homogeneous and isotropic at large scales. Section 4 devises a simple strategy to solve the 
coincidence problem. Finally, Section 5 summarizes and discusses our results. 

We focus our attention on spatially flat Friedmann-Robertson- Walker models. As usual, 
a subscript zero means that the corresponding quantity is to be evaluated at present time. 



II. THREE CHOICES FOR THE INTERACTION TERM 

Because of the unknown nature of both DM and DE it does not seem feasible nowadays to 
derive an expression for Q on first principles. The alternative is to propose phenomenological 
expressions based on arguments of simplicity and reasonableness. Firstly, note that Q must 
be small and positive -otherwise the DE could not dominate the expansion today. On the 
other hand, if it were large and negative, DE would have dominated the expansion practically 
from the outset and galaxies couldn't have condensed. 

Looking at the left hand side of Eqs. ([T]) it becomes apparent that Q must be a function 
of the energy densities multiplied by a quantity with units of inverse of time. For the latter 
the obvious choice is the Hubble factor, so Q = Q{Hpm, Hp^). By power law expanding 
this function and keeping just the first term we get Q ~ X^H + X^H p^. In view of 
the lack of information about the coupling it is advisable to retain just one parameter (the 
less observational parameters when fitting to observational data, the better). Then, three 



choices follow, namely: Xm = Xx, Xx = 0, and Am = 0. Other choices, as Q oc pmPx 



are 
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also possible but we will not deal with them here. 

We will therefore consider the following three expressions for the coupling, 



Q = 3aH{pm + Pct>) , Q = ^PH pra , and Q = 3r]Hp^ , (4) 

where the phenomenological parameters, a, j3, and 77, are dimensionless, positive-definite 
small constants. Henceforward, we will call the corresponding models alpha, beta, and eta, 
respectively. 

Model alpha fits very well data from SN la, CMB, and large scale structure formation 



provided a < 2.3 x 10 



21 
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251]. Further, this model possesses the remarkable property 
that the ratio r tends to a stationary but unstable value at early times, r^, and to a 
stationary and stable value, r~, (an attractor) at late times. This follows from inserting the 
corresponding expression for Q in the right hand side of Eq. ([2]) and setting r to zero. For 
w = constant the stationary solutions of the resulting quadratic equation are: 



r 



± 



-l + 2b±2Jb{b-l) , b =-—>!. (5) 

4a 

As it can be checked by slightly perturbing the right hand side of Eq.([2]), the stationary 



solution rt results unstable while is stable 



26 



27| . The general solution of Eq. 



(6) 



interpolates between r+ and r~ . Here, x = (a/a*)~^, with p = 12 a^Jb{b — 1), and 
denotes the scale factor at which r takes the arithmetic medium value (r+ + r^)/2. In the 
range r~ < r < the function r{x) decreases monotonously. Consequently, as the Universe 
expands, r{x) smoothly evolves from r+ to the attractor solution rj. The transition from 
one asymptotic solution to the other occurred only recently -see left panels of Figs. [H and 
[2]- so one can take r ~ during a fairly large part of the history of the Universe. Finally, 
the constraint ~ constant implies that r and a are not independent, but verify the 
relationship a{rf + 1)^ =| w \ r+, so the product arf w \ is of order unity. 

As indicated above, this choice of Q yields a constant but unstable ratio at early times 
and constant and stable lower ratio at late times -left panels of Figs. [1] and [21 It is hard to 
imagine a simpler expression for the coupling entailing these two properties. Likewise, the 
aforesaid expression can be re-interpreted as implying an effective exponential potential for 



the dark energy field at late times and a power-law at early times 



24 



26| 
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FIG. 1: From left to right, evolution of the ratio r = pm/ p<f) with redshift for models alpha, beta, 
and eta. For all of them r either tends to a constant or varies very slowly at small redshift. In 
drawing the curves we have fixed tq = 3/7 and w = —0.9. 
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FIG. 2: Same as Fig.l except for that in this case w = —1.1. 

Near r r~ the balance Eqs. ([1]) can be approximated by 

1 dp,^ _ 1 - a(l + l/r-) 
pm dr ~ a{r+ - r^){r - r-) ' 
1 dp^ l + w + a(l + r^ 



P0 dr a{r+ - r-){r - r-) ' 
For w ~ constant, these equations can be integrated to 

Pm oc a-4i-"(i+iA'^')] , p^ oc a-3[i+«,+.(i+r. 



(7) 



(8) 



Notice that the condition r\r=r- = implies that the exponents in the energy densities, 
Eq. ([8]), coincide. Interestingly, these results are not only valid when the dark energy is 
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a quintessence field (i.e., —l<w< —1/3), they also apply when the dark energy is of 
phantom type (i.e., w < —1), either a scalar field with the "wrong sign" for the kinetic 
energy term, a ^-essence field, or a tachyon field 



In the vicinity of the attractor, dark matter and dark energy dominate the expansion 
whence the scale factor can be written as 

^ ^ ^(2/3)[l+u>+a(l+rr)]-l_ 

It is easy to show that, at early times, + {wa~^ + 2)r + 1 = and^ the constant ratio 
is r+ = r{z oo) ~ a^^ while, at late times, r~ = r{z — > — 1) ~ a [27j]. These results 
mean a conceptual advance with respect to the ACDM model. If the present acceleration 
of the Universe is generated by a cosmo logical constant, then the initial values of dark 
matter and dark energy have to be fine-tuned by 96 orders of magnitude at Planck time, 
i.e., ritpianck) = ^m/^A ~ 0(10^^). In the model of above this ratio is fixed by a. Thus, 
the initial condition problem is significantly alleviated. For example, if a ~ 10"'^, the initial 
dark matter to dark energy ratio comes to be r(tpianck) ~ 0(10"^). 

Model beta does not share the nice feature of a late attractor. However, in this case, for 
constant w one has {29] 

'^^/5ro-(l + ^)-3«[/3(l + ro)+H ' 

where ^ = —w — j3 > 0. Hence, the evolution of r is slower than in the ACDM model. 
There is a monotonous decrease of r from values higher than unity at early times to a 
nearly constant, value at present time (fixed as 3/7) irrespective of whether dark energy is a 
quintessence field (middle panel of Fig. [T]) or a phantom field (middle panel of Fig. [2]). For 
vanishing scale factor r tends to the finite, constant value ^//3 (however, the model should 
not be extrapolated to such an early stage). Thus, in spite of the absence of a late time 
attractor this choice of Q also significantly alleviates the coincidence problem as the energy 
density ratio varies very slowly at present time, i.e., | (r/r)o |< Hq. This feature is shared 
by model eta -see right panel in Figs. [T] and [H 

Altogether, all three models considerably alleviate the coincidence problem but there is 
no guarantee that any of them solves it in full. A successful model should also predict 
tq ~ C>{1). To the best of our knowledge, no model in the market does this. We believe 
that, for the time being, the present value of r must be considered an input parameter much 



in the same way as the current value of the Hubble function, Hq, the temperature of the 
CMB, the age of the universe, or the ratio between the number of baryons and photons, are 
-only that the four latter are not laden with any puzzle. 



III. COMPARISON WITH OBSERVATIONAL DATA 

We wish to explore whether the interacting, models alpha, beta, and eta, are consistent 
with the results of Daly et al. j3] • These are based in the determination of the dimensionless 
coordinate distance, y{z) = Hq (aQr), and their two first derivatives with respect to redshift. 



of supernovae type la (SN la) and radio-galaxies (RG) in t 
The supernovae samples contains 182 SN la of Riess et al. 



le redshift interval < 2; < 1.8. 



30|, 192 of Davis et al. (31| and 



115 of Astier et al. 321] though there is some overlap between the different samples. The 
RG sample comprehend 19 of Daly et al. [33] plus 11 additional RG. 

It is to be stressed that the determination of the dimensionless coordinate distances 
does not involve any cosmological model, just the Friedmann-Robertson- Walker (FRW) 
metric. The first and second derivatives of the dimensionless coordinate distance with respect 
to redshift can be used to build a model-independent determination of the dimensionless 
expansion rate, E{z) = H{z)/Hq, and the deceleration parameter, q{z) = —a/{aH'^), see 



34]. 



Indeed, in virtue of the relation between redshift z and the scale factor 1 + z = ao/a{t), the 
definition of the Hubble rate, H = a/a, and the relationship dt = —a{t) dr, valid in spatially 
fiat FRW spacetimes for photons fiying from their source to the observer, one follows: 

. ^ H(z) 1 , , 

where a prime indicates derivation with respect to z. From last expression the data y(z) 
can be used to empirically determine y\z^ and E{^z\ Differentiating Eq. fllOl) with respect 
to time and using the definition of g(z) leads to 

g(z) = -l-(l + z)^. (11) 

Likewise, differentiating Eq.Q with respect to time and using Eqs. ffTOl) and ffTTl) produces 



P^{z) = P^{z = 0) 



E'^iz) \ r 2q{z) - 1 



E'^{z = 0) 



2q{z = 0) 



(12) 
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Thus, the dimensionless quantities E{z), q{z) and P^{z) can be constructed directly from 
y' and y" js^]. 

In virtue of Eqs.([ID,([2]), ^ and ^ we get 

y'i^) = T7? , (13) 



and 



2{l + z) 



(14) 



2 

These expressions are vahd for all three models of above. Here f2mo = ^bh^ ^^"^ ^</'0 ~ 

2 

"^^^2° denote the present day values of the density parameters of DM and DE, respectively. 
Numerical integration of Eq. (fT3ll yields the dimensionless coordinate distance, y{z) = 
y\z) dz. In doing this, w = constant is assumed. 

Next, we consider the three expressions for the source term Q introduced above. 

(i) For model alpha (interaction term Q = 3a H [p^ + Pm))-, we find that 

^i{z) = {l + zY\ and = (1 + 2)"+"^ , (15) 

with 77, = I [2 + w — ^Jw{w + 4q;)] , and m = 3 \Jw{w + 4a). 

(ii) For model beta {Q = SP H pm) we obtain 

3i(;.) = (1 + zr , and ^32 = (1 + - --^ ^ 53i(z), (16) 

with n = 3(1 -/?). 

(uz) Likewise, for model eta {Q = Sr] H p^) 

3i(z) = (l + zf -^-^532(^) and ^2 = {l + zf'^'+-+^\ (17) 

In the absence of interaction expressions ffTB]l . ffTB]) . and ffTTl) reduce to Eqs. (1) and (2) 



of Daly et al. 



'221. 



The evolution of the quantities y, y' and y" in terms of z are depicted in Figs. [3] through 
m for the three models. We have used different values of the parameters w, a, /3, and rj, 
and assumed the standard values VL^q = 0.3 and VL^q = 0.7. Also, the said evolutions of y, 
y' and y" are compared with the flat ACDM model with the quoted values of VL^o and VL^q. 

The data and analysis from the values of y{z), y'{z) and y"{z) are used to obtain E{z), 
q{z) and P^{z) / P^{z = 0). The results are shown in Figs. [6] through [8] and, again, contrasted 
with the standard ACDM model. 



9 



Comparison of Figs. [3] through [7] with the corresponding figures of Ref. [22j reveals 
that the three models, alpha, beta, and eta, are consistent with the analysis of Daly et 
al. which is based on different samples of SN la and 30 RG. It is worthy of mention that, 
except for adopting the Friedmann-Robertson- Walker metric, the study is independent of 
any cosmological assumption. We, therefore, conclude this section by stressing that the 
three models are consistent with the known properties of dark energy. 

IV. A STRATEGY TO SOLVE THE COINCIDENCE PROBLEM 

As mentioned above, a cosmological model that predicts a constant or slowly evolving 
ratio r at late times certainly alleviates the coincidence problem but this is not enough. The 
model also must predict that tq ~ 0{1). Our strategy is very simple; it rests on the natural 
assumption that r{z) and H{z) must be closely interrelated. Therefore, if one express r in 
terms of H the problem of explaining why tq is of order unity reduces to the less pressing 
problem of explaining the present value of H . This is interesting because while we are 
not aware of model-independent data of r at different redshifts, two sets of observational 
data values -though scarce- of H at various redshifts are now available. One data set was 
compiled by computing differential ages of passively evolving galaxies in the redshift range 
0.1 < 2; < 1.8 [35|. The other data set is based on measurements of supernovae type la (SN 
la) and 30 radio galaxies up to redshift 1.2 {22! . 

To express r as a function of H we use the identity f = H dr/dH, with 

H = -'^{p^ + p, + P,) = -l'-±^H\ (18) 



in Eq. Then, the latter becomes 

dr Q 
dH^H 

where 



l + r 

3 = -2r — 



1 + w + r 



«:2Q(r + l)2 



(19) 



(20) 



Because of we expect dr/dH > (a very reasonable assumption since H(t < to) > Hq as 
observation tells us) the inequality 
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is to be satisfied. 

We note in passing that for tlie particular case of tlie ACDM model (i.e., w = —1 and 
g = 0) one has H = Ho [(1 + r)/(l + tq)]^/^. 

A. Model alpha 



For this model, Eq. ( l20l) reduces to 



$5 = -2r 



1 +r 



1 + w + r 



w + a 



[r + 1)^ 



(22) 



where we used Eq. (|3]). The condition given by Eq. fl2Tl) leads to restriction li^l > a (1 + 
r)^/r. 

Using last equation in ( fT9l) and integrating, the Hubble function can be expressed in 
terms of r as 



i7(r) = ifoexp[/(r)-J(ro)], 



(23) 



where /(r) stands for the real part of J(r) with 



/(r) = - In(H-r) — - ln[— w r — a(l + r)^] 



1 



w + 2 



tan 



2a(l + r) + tt; 
a/— w (4a w) 



• (24) 



2 ■>/— w(4a + w) 

Figure [9] shows the dependence of the Hubble function on the densities ratio for a = 
10~^ and two values of the equation of state parameter, w. Likewise, we have plotted the 
prediction of the ACDM model. In all the cases vq = 3/7. For other a values compatible 
with the restriction a < 2.3 x 10^^ set by the Wilkinson Microwave Anisotropy Probe 3 yr 
experiment the corresponding graphs (not shown) are rather similar. 

Unfortunately, as said above, no model-independent data of r at different values of H 
(or z) exists whereby, for the time being, we cannot directly contrast this prediction with 
observation. This is why we turn to determine the dependence of the Hubble factor with 
redshift and compare it with the two available observational data sets of H vs z. The 
first model-independent data set was obtained from the study of the differential ages of 32 
-carefully selected- passively evolving galaxies in the redshift range 0.1 < z < 1.8 jssj. 
The age of each galaxy was found by constraining the age of its older stars with the use of 
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synthetic stellar population models. The differential ages roughly yields dzjdt, then H(z) 
is given hj H = —dz/[{l + z) dt] -see right panel of Fig. 1 in Ref. 35j. The second 
model-independent data set (see lower panels in Fig. 9 of Daly et al. {22]) was obtained by 
applying the model-independent analysis method of Ref. 34 |to the coordinate distances of 
192 SN la of Davis et al. 3l| and 30 radio galaxies of Ref. 33|. 

To obtain the expression for H(z) of this model we combine Eq. (123!) with the integral of 
(12]) in terms of z which is 



r{z) = Re 



1 

2a 



-w 



7 tan ( A — - 7 ln(l + z) 



1 



(25) 



where Re specifies the real part of the quantity in curly parenthesis, A = tan ^ [{w + 2a(l + 



ro))/7] , and 7 = ^/-w {w + 4 a;). 

Figures [TO] and [TT] show the dependence of the Hubble expansion rate on redshift predicted 
by the model for a = 10~^ and two values of the equation of state parameter, w. For 
comparison, we have also plotted the prediction of the ACDM model. The data in Fig. [10] 
are borrowed from Simon et al. 35|] (full circles) and Daly et al. (full diamonds), and in 
Fig. [Ufrom Daly et al. (cf. lower panels of figure 9 in Ref. 22|). As can be seen in both 
figures, the fit to data does not vary significantly between graphs. Given the scarcity of data 
it is not worthwhile to compare the of the different curves. 

Figures [9] [10], and [11] taken together suggest that this model may be a solution to the 
coincidence problem. It fits the available data sets of the evolution of the Hubble factor 
not less well than the ACDM model. Unfortunately, due to the large la errors (especially 
for z > 0.5) they cannot discriminate between these two models. A like situation occurs 
for the beta and eta models detailed below. In consequence, more abundant and accurate 
observational information about the aforesaid evolution is needed in order decide whether 
any of the three models can really solve the coincidence problem. 



B. Model beta 



This model has been considered in 
reduces to 



36| and (37|. Equation ([20]), with the help of 



c> = _2r _ [w + (3 (r + 1)] 



1 + w + r 

Now, the condition set by Eq. ([2T]) boils down to \w\ > P {1 + r). 



(26) 
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In virtue of Eqs. (fT9l) and (126!) the Hubble function can be cast as 

ff(r) = i/„^. (27) 

where 

/(r) = v^(r+7) r^^-'-m^iP+u,)] ^^(^ ^ ^) ^ ^p-;3)/[2(/3+«,)] _ (28) 

Numerically, the dependence of the Hubble expansion rate upon r is very close to that 
of model alpha for (3 values similar to a. This is why we do not show it here. 

The dependence of the Hubble expansion rate on redshift can be obtained by combining 
Eq. (!27j) with ([9]). It is found that the corresponding graphs for H{z) are rather similar to 
those in Figs. [10] and [11] whereby we do not depict them here. 

C. Model eta 

This model has been studied in I23I. In this instance, Eq. fl2U]) becomes 

(r + 1)' 



1 + r 

3 = -2r 



W + 7] 



(29) 



1 + w + r 

and the condition set by Eq. f[2T]) simplifies to \w\ r > ri{l + r). With the help of Eqs.ffW 
and ( 129]) the Hubble function, H{r), is given formally by ( 1271) but this time 



I{r) = y'{l + r) [r] + r {w + rj)] '^(^+^) . (30) 
By combining the expression for H{r) with the integral of ([2|) in terms of z, 

^(^) = f \ ^ + v)+v{l-{l + zf-^^>^)} , (31) 

[w + ri) ^ 

one follows H{z). Again, the graphs of H{z) are very close to those of Figs. [TO] and [TT] and 
we do not show them here. 



V. DISCUSSION 

Nowadays the leading candidate for dark energy is the popular ACDM model. It is 
simple, it has just one free parameter -the vacuum energy density- and fits reasonably well 
most (if not all) observational data. However, it presents the puzzle (apparent or not) that 
we happen to live at a very special epoch: the transient epoch at which the ratio between 
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the densities of matter and vacuum energy is of order unity. To solve this one must look 
for alternative models in which the ratio r stays constant or varies very slowly, around the 
present time, with respect to the Universe expansion. Interacting models (with Q > 0, see 
Eq- (ED) possess this property. But it does not suffice. The successful model must also 
predict tq ~ 

In this paper we considered three interacting models that present the ffist characteristic 
and may also have the second one. To decide on the latter one must contrast the r{z) 
function predicted by each model with observation. Regrettably, as far as we know, model- 
independent data about r{z) do not exist. However, albeit with wide uncertainties, there 



exist model-independent data of the evolution of the Hubble factor up to z = 1.8 [22|, |35 |. 
Therefore, we have numerically determined H{r) and H{z) for each model and compared 
the latter with these observational data. All three models show compatibility with the data 
but only in a degree similar to the ACDM model. Future, more accurate and extensive 
measurements will likely tell us whether the latter model fits the data better or worse than 
interacting models (the ones presented in this paper or maybe others). If, in the end, the 
ACDM model makes the best fit, then we may conclude that the coincidence problem is not 
a problem any more; it is simply a coincidence and therefore not a problem. 
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FIG. 3: From top to bottom, dimensionless coordinate distance, y{z) = Hq {aof), in terms of the 
redshift for models alpha, beta, and eta. For comparison, the prediction of the ACDM model is 
also shown. 17 




FIG. 4: From top to bottom, first derivative of the dimensionless coordinate distance in terms of 
redshift for models alpha, beta, and eta. For comparison, the prediction of the ACDM model is 
also shown. 18 
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FIG. 5: From top to bottom, second derivative of the dimensionless coordinate distance in terms 
of redshift for models alpha, beta, and eta. For comparison, the prediction of the ACDM model is 
also shown. 
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FIG. 6: Prom top to bottom, evolution of the Hubble factor for models alpha, beta, and eta. For 
comparison, the prediction of the ACDM model is also shown. 
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FIG. 7: Prom top to bottom, evolution of the deceleration parameter for models alpha, beta, and 

eta. For comparison, the prediction of the ACDM model is also shown. 
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FIG. 8: From left to right, evolution of the ratio P^{z) / Pfj){z = 0) with redshift for models alph 
beta, and eta. 
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FIG. 9: Evolution H vs the ratio r = Pm/P4> as given by Eq. (p3]) -model alpha- with a = 10~^ 
for w = —0.9 (quintessence, top line) and w = —1.1 (phantom, bottom line). Also shown is the 
prediction of the ACDM model (middle line). In drawing the curves we have set rg = 3/7. 
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FIG. 10: Evolution H vs z with a = 10~^ for w = —0.9 (quintessence) and w = —1.1 (phantom). 
Also shown is the prediction of the ACDM model (solid line). In all the cases we have fixed 
ro = 3/7 and Hq = 71km/s/Mpc. The data points with their la error bars are borrowed from 



Simon et al. Ref. |^5j (full circles) and table 2 of Ref. 22] (full diamonds) 
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?1G. 11: Same as Fig. [TOl The upper panel shows the combined sample of 132 SN la of Ref. 



31[ and 30 radio galaxies of Ref. 



33( 1 ■ The bottom panel shows the 30 radio galaxies only. The 



data points with their la error intervals and the best-fit curve (big solid line) in both panels are 
borrowed from Daly et al., Ref. [22]. 
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